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We theoretically investigate the confinement-induced resonance for quasi-one-dimensional quan- 
tum systems under transversely anisotropic confinement, using a two-body s-wave scattering model 
in the zero-energy collision limit. We predict a single resonance for any transverse anisotropy, whose 
position shows a slight downshift with increasing anisotropy. We compare our prediction with the 
recent experimental result by Haller et al. [Phys. Rev. Lett. 104, 153203 (2010)], in which two 
resonances are observed in the presence of transverse anisotropy. The discrepancy between theory 
and experiment remains to be resolved. 



C/2 

a- 
a ■ 



o 



(N 
> 

(N 

o 
o 



X 



I. INTRODUCTION 



Recently ultracold low-dimensional atomic gases have 
attracted a great deal of interest as they show unique 
quantum signatures not encountered in three dimen- 
sion (3D) [l|. For example, a one- dimensional (ID) 
gas of bosonic atoms with strongly repulsive interpar- 
ticle interactions acquires fermionic properties [6]. A 
ID gas of fermionic atoms may exhibit exotic inhomoge- 
neous superfluidity Q-Sl and cluster- pairing Q- In two 
dimensions, the Berezinskii-Kosterlitz-Thouless (BKT) 
phase transition is very different from the usual finite- 
temperature phase transitions [7]. Experimentally, ul- 
tracold low-dimensional atomic systems are conveniently 
accessible upon introducing tight confinement via optical 
lattices that remove one or two spatial degrees of free- 
dom [l|. Furthermore, the interactions between atoms 
can also be tuned precisely and arbitrarily by a colli- 
sional Feshbach resonance jaQ- Using these means, low- 
dimensional atomic systems now provide an ideal exper- 
imental model to test fundamental many-body physics, 
and are much simpler than corresponding systems in con- 
densed matter physics. 

Confinement-induced resonance (CIR) is an intriguing 
phenomenon found when one relates the effective ID cou- 
pling constant of quasi-lD quantum gases, gio, to the 
3D s-wave scattering length, a^^o |10l - [la |. Considering a 
tight transverse confinement (o;^ ~ '^x ~ '^y) with the 
characteristic length scale a±^ — yjh/{muj±^, Olshanii 
showed that the effective ID coupling strength diverges at 
a particular ratio of a^D/a^ [10]. The underlying physics 
of the CIR is very similar to the well-known Feshbach res- 
onance [11], if we assume that under tight confinement 
the ground-state transverse mode and the other trans- 
verse modes play respectively the roles of the scatter- 
ing "open" channel and "closed" channels. The scattering 
state of two colliding atoms in the ground-state trans- 
verse mode can be brought into resonance with a bound 
molecular state in other high-lying transverse modes by 
tuning the ratio of aso/o-^- This causes a confinement 



induced resonance. 

This type of novel resonance was recently confirmed in 
several ultracold atom experiments [16 20]. In particu- 
lar, the properties of CIR for an ultracold quantum gas of 
Cs atoms were studied systematically by measuring atom 
loss and heating rate under transversely anisotropic con- 
finement near a rather narrow Feshbach resonance [2C|. 
Upon increasing the transverse anisotropy, two CIRs were 
observed, splitting from the single CIR in the isotropic 
limit. 

In this paper, we theoretically investigate the CIR un- 
der transversely anisotropic confinement, using a zero- 
energy pseudopotential approach [10], which is a useful 
generalization of Olshanii's method. We find a single res- 
onance in the scattering amplitude, except with a down- 
shift of the resonance position with increasing transverse 
anisotropy. We calculate also the bound state of the de- 
coupled excited transverse manifold and find only one 
such bound state exists. Taken together, we believe that 
there is only one CIR in an anisotropic two-body s-wave 
scattering model at zero-energy limit. The two (split- 
ting) resonances found experimentally therefore must be 
caused by some other reasons beyond the simple s-wave, 
zero-energy scattering model. 

In the following, we first present the essence of s-wave 
scattering approach and calculate the bound state of the 
excited transverse manifold (Sec. II) and then discuss in 
detail the CIR as a function of the transverse anisotropy 
(Sec. III). Sec. IV is devoted to conclusions and further 
remarks. 



II. CIR UNDER TRANSVERSELY 
ANISOTROPIC CONFINEMENT 

A. The scattering amplitude 

In this section, we calculate the effective ID coupling 
constant under transversely anisotropic confinement by 
generalizing the zero-e nerg y s-wave scattering approach 
developed by Olshanii [1(||. Let us consider two atoms 



in a harmonic trap with weak confinement in z-axis and 
tight confinement in the transverse direction such that 
LOx = Tj^y 3> ujz- The transverse anisotropy is charac- 
terized by the ratio 77 — uix/ujy ■ The motion along the 
z-axis is approximately free. To consider the scattering in 
this direction, we therefore set ujz — 0. Owing to the sep- 
arability of center-of-mass motion and relative motion in 
harmonic traps, the low-energy s-scattering along z-axis 
is described by the following single-channel Hamiltonian: 






(1) 



where 
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Here, gac = iirh'^asD/'m is the 3D coupling constant, 
and /i = m/2 is the reduced mass for the relative mo- 
tion. The effective ID coupling constant giD is obtained 
by matching the scattering amplitude of the exact 3D 
solution ^(r) of Eq. ([l} to that of a ID Hamiltonian 






(3) 



The latter scattering amplitude feikz), with low scatter- 
ing energy E = ?i^fc^/2/i, is then given by, 

^{x = 0, y = 0, z) oc e*'^^^ + /e(fc.)e^'=^l^l (4) 

for |z| — >■ 00 , and the effective ID coupling constant guj 
can be determined by using 
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Next, we shall calculate the 3D exact solution ^(r), 
from which we can extract the effective ID scattering 
amplitude /e(fcz) at the collision energy E = ti^k'^/2^ + 
huJx/2 + fkjjy/2, when kz — >■ 0. Using the complete set of 
eigenstates of the non-interacting Hamiltonian, the (un- 
normalized) scattering wavefunction can be expanded as. 
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(6) 

where (/)„(0 = {y/^2''n\)-^/^e-^^^^Hn{0 is the n-th 
eigenstate of a ID harmonic oscillator (n = 0,2,4,...), 
and d = ^Jhj (utOy) is the oscillator length in the j/-axis. 
In the sumination, the integers ni,n2 can be restricted 
to even numbers, as s-wave scattering preserves parity. 
The expansion coefficients A„j„2(^) may be determined 
by substituting the wavefunction ^ into the Schrodinger 
equation Hrei'^i'^) — E'^ir) and then projecting both 
sides of the equation onto noninteracting states (i.e., the 
expansion basis) [21-23]. Owing to the contact interac- 
tion, the expansion coefficient takes the form. 
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depending on a single parameter 
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The summation over ni, 77,2 and k in the scattering wave- 
function ^(r) can then be performed [22,[23|, yielding 
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where the function A[^, e] takes the form 
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The parameter J- is to be determined and it should satisfy 
the boundary condition, Eq. ([S]). We thus have to check 
the asymptotic behavior of the wavefunction (]§]) as z — >■ 
0+, or to check A[^, e] as ^ — >■ 0. After some algebra, one 
finds that. 
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At low energy, the zero-order term of £ (e) has the form 
C (e) = — C + £ (e), where the constant C is defined as 
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Using the above asymptotic form of A[f , e], the parameter 
T can be easily determined, and hence the effective ID 
scattering amplitude by comparing Eq ® and Eq. Q. 
We finally arrive at, 
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In the limit of an isotropic transverse confinement (77 
1), the constant C is 
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The function C (e) has the form 
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Using the integral form of the Riemann zeta function 
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we find that 
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Eqs. (fT5|) and (fT5)) reproduce the well-known result by 
Olshanii [10]. 

We can easily obtain the ID coupling constant gin 
from Eq.© and Hi 
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which diverges at {a^^J/ay) = y/2/C , where ay = 
^hjraujy — d/\/2 . Therefore, there is only one CIR, 
whose position depends on the ratio 77 — ujx/ujy for a 
given length scale ay. At 77 = 1, the constant in Eq. P^ 
can be calculated analytically, giving C,,=i = — C(l/2) — 
1.4603.... This recovers the well-known result given by 
Olshanii 01 . 
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Figure 1: (Color online) The bound states of the total relative 
Hamiltonian Hrei (the solid curve) and the excited Hamilto- 
nian He (the dotted curve) with 77 = 1 . 
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Figure 2: (Color online) The same as Fig[l]but for rj = 2 



B. Bound state of the excited transverse manifold 

Bergeman, et al. [ll| show that the physical origin of 
the CIR is a zero-energy Feshbach resonance under tight 
transverse confinement. The ground transverse mode 
and the manifold of excited modes play the roles of the 
open channel and the closed channel, respectively. If 
a bound state exists in the closed channel, a resonance 
will occur when the bound-state energy is equal to the 
incident energy of two colliding atoms in the ground 
transverse mode of the total Hamiltonian. In this sec- 
tion, we calculate the bound state of the decoupled ex- 
cited manifold. We find that a resonace will appear 
at a\jjjay = v2/C, consistent with the calculations in 
Sec|^ 

Using a similar method to Sec lII Al the wavefunc- 
tion of the total relative Hamiltonian Hj-ei can be ex- 
panded in terms of the eigenstates of the non-interacting 



Hamiltonian, Eq.®. After substituting Eq.([6]) into the 
Schroedinger equation and using the normalization and 
orthogonality of the wavefunctions, we arrive at. 
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where 



H(f) 



dkz 



00 00 

EE 



Jk^zi* 



K, (0) C. (0) K^ {x) K, (y) 



r^ 27r 

-°° ni=0n2=0 



\d'^k1 + nxr\ + n2 — e 

(21) 
and e — E/fkjy ~ r\l1 — 1/2 . Using 
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and 
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Eq. (|2ip can be written as 
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where Pg = \g) {g\ is the projection operator of the trans- 
verse ground mode \g) . We can also calculate the bound- 
state energy of the excited Hamiltonian He by removing 
the Til = 712 = term from the summation of the func- 
tion Hif) , i.e., Eg. ipT]) . Then the bound-state energy 
of the excited Hamiltonian Hg should be determined by, 
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for e < . As r ^^ , only the small t terms dominate 
the integral, so that 
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for e < . Consequently, we can redefine the regular part 
of He {r) as 



V' 



2Trd^r ■ 
Hence, the function H (r) can be rewritten as. 
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and the regular part J- (e, f) has the form, 
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Then, using Ea. ([5T|) and ([55)) . we arrive at 
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Next, from Eq.(|20p. we obtain. 



t3/2 
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This determines the bound-state energy of the excited 
Hamiltonian He ■ 

For isotropic confinement {ij — 1), it is obvious that. 
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for e < . Eq. ipSj) and (^5]) determine the bound-state 
energy of the total relative Hamiltonian depending on the 
interaction strength, which is denoted by the 3D scatter- 
ing length a^D ■ 

If the transverse ground state is projected out from the 
Hilbert space, we can define the excited Hamiltonian as. 



He = Hr 



rgTl^.el^g 



(30) 



(35) 



Hence, from Eg. ([25)) and ([M]) . the bound state of the 
excited Hamiltonian He is 2f)w higher than that of the 
total relative Hamiltonian Hrei , which is consistent with 
the conclusion of Bergeman |ll| . 

The bound states of the total relative Hamiltonian Hrei 
and the excited Hamiltonian He are plotted as functions 
of the 3D scattering length a^o in Fig[T] and Fig[2] with 
transverse aspect ratios ?/ = 1 and rj = 2 for isotropic 
and the anisotropic confinement, respectively. It is ob- 
vous that there is only one bound state for the excited 
Hamiltonian He both for isotropic and anisotropic con- 
finement. A confinement-induced resonance occurs when 
this bound state becomes degenerate with the thresh- 
old of the transverse ground mode of the total relative 
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Figure 3: (Color online) The effective ID coupling constant, 
giD (in units of f? /{may)), as a function of a^o /ay at different 
transverse anisotropy 77 = 1 (solid line), 1.5 (dashed line), and 
3 (dot-dashed line). 



Hamiltoiiian. The position of the CIR in terms of the 3D 
scattering length is determined by Ea. ([M)) with e = 0, 
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This result is consistent with that in the Sec lII Al which 

is (45 /a,) = V2/C. 

III. RESULTS AND DISCUSSIONS 

We plot in Fig[3]the effective ID coupling constant as a 
function of the 3D s-wave scattering length at 77 = 1, 1.5, 
and 3. In this parameter range, with increasing trans- 
verse anisotropy, the resonance position of CIR shifts 
from (a3D/aj^)[^^\ ^ \/2C~ii = 0.9684... to lower val- 
ues. For the experimental situation with a 3D scattering 
length tuned by a magnetic Feshbach resonance, this de- 
crease corresponds to a shift to the deep BEC side of the 
Feshbach resonance. 

In a wider parameter space, however, we observe an 
interesting non-monotonic behavior of the resonance po- 
sition as a function of the ratio 77. As shown in FigUl 
at sufficiently large (or small) values of the transverse 

anisotropy 77, the resonance position a^jj increases with 
increasing (or decreasing) transverse anisotropy, and di- 
verges at a critical frequency ratio 770 — 18.551 (or 

— 1 ( R) 

rjQ ~ 0.0539). For Og^, , this corresponds to a bend- 
back to the magnetic Feshbach resonance from the BEC 
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Figure 4: (Color online) Resonance positions of CIR, a]^jj /ay, 
as a function of the transverse anisotropy 77 — cux/ij^y For 
sufficient large or small values of transverse anisotropy, the 
quasi-lD system crosses over to the quasi-2D regime. 



limit. On further increasing (or decreasing) the trans- 
verse anisotropy, the quasi-lD quantum gas crosses over 
to a pancake-shaped, horizontally-oriented 2D system, 
which would lead to a negative value for a^J ■ This ob- 
servation is in agreement with theoretical calculations [1J| 
showing that for a quasi-2D system, the s-wave CIR oc- 
curs always at a negative 3D scattering length {aso < 0) 
or on the BCS side of the magnetic Feshbach resonance. 
We now compare our theoretical results with the ex- 
perimental data [20j. This experiment uses a narrow 
Feshbach resonance in ^^^Cs ^], together with a nona- 
diabatic jump to a particular magnetic field to initiate 
the resonance with a finite kinetic energy. As a result, 
one may expect some differences to the simplest single- 
channel model used here. The single-channel model 
does not include the molecular field corresponding to 
the bound molecule state. It is known from both mean- 
field [231 and beyond mean- field [2a| calculations that, 
the single-channel model describes only a broad Fesh- 
bach resonance. As illustrated in FiglSJ our prediction 
of a^J is just between the two CIRs observed in the ex- 
periment, and disagrees with either of them. In fact, the 
observation of a single CIR s-wave resonance in our the- 
oretical calculation can be easily understood. Recall that 
the CIR is a kind of Feshbach resonance, which occurs 
when the scattering state in an open channel becomes 
degenerate in energy with the bound states in the closed 
channel. The number of CIRs is therefore simply equal 
to the number of bound states in the closed channel. For 
a quasi-lD quantum gas, the excited transverse modes 
are grouped together and regarded as closed channels 
(see, for example, ref. [11| for a detailed discussion of 
closed channels). In our model, the sub-Hilbert space of 
the excited transverse modes can only sustain a single 
bound state, no matter how large large the transverse 
anisotropy is. The experimental observation of two res- 



onances must therefore have a different origin than the 
simple zero-energy s-wave scattering picture pursued in 
this work. 

Interestingly, we obtain a good agreement with exper- 
imental data if we include in the resonance condition 
an additional factor of y/rj. As shown in Fig. 5, the 
CIRi and CIR2 resonances are well reproduced theoret- 



ically if we use the conditions, 



ho 






/ay 



2ri/C and 



/ay 



/2/77/C, respectively. This additional factor 



however can not justified. 

It should be note that, by keeping ujx invariant and 
weakening the confinement ujy - which therefore gives an 
alternative way to increase the value of 77 - the observed 
CIRi displays additional structure. It further splits into 
many sub- resonances at r/ > 1.5 and finally disappears 
at sufficiently large values of transverse anisotropy. At 
the same time the CIR2 persists. Also, in the limit of a 
2D system (77 — > cxo), the resonance position of the sur- 
viving CIR2 is at a positive 3D scattering length. These 
experimental observations also disagree with our simple 
zero-energy s-wave scattering calculation, and with pre- 
vious calculations^] of 2D s-wave CIR. One possibility 
is that many-body physics of the quantum bosonic gas 
may come into play. A likely explanation therefore is 
that the CIRs observed in this experiment are associ- 
ated with the failure of Olshanii's model [10| which is 
only a zero-energy single-channel s-wave description of 
the inter-atomic interactions, and the experimental con- 
ditions are different to this simple model. An advanced 
theory of CIRs need to be introduced, e.g., multi-channel 
scattering theory including molecular channels, non-zero 
collision energy theory, or many-body effects. 



IV. CONCLUSIONS AND REMARKS 

In conclusion, we have presented a theoretical cal- 
culation of confinement-induced resonance (CIR) under 
transversely anisotropic confinement, by extending the 
standard zero-energy s-wave scattering approach to in- 
corporate transverse anisotropy. We have theoretically 
found a single CIR and have explained its physical ori- 
gin. The position of our predicted CIR disaggrees with 
either of the two CIRs measured most recently in a 
quasi-lD quantum bosonic gas of ^^^Cs atoms. There- 
fore, the CIRs observed in the experiment cannot be 
explained simply by the zero-energy s-wave scattering 
theory. There are also several discrepancies when the 
quasi-lD system crosses over to the quasi-2D regime, 
which may arise either from the many-body physics of 



Cs atoms, or the residual molecular fraction and angular 
momentum of the magnetic Feshbach resonance. To take 
these aspects into account, we look forward to improving 
the current s-wave scattering calculation by replacing the 
two-body scattering matrix with its many-body counter- 
part and/or adopting a two-channel model to describe 
the magnetic Feshbach resonance. 
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Figure 5: (Color online) Resonance positions of CIRs, ob- 
tained by solving a\,D /°'-v ~ V^/C, as a function of the 
transverse anisotropy 77 = uj^/^y The experimental results 
(symbols) are compared with our zero-energy s-wave scat- 
tering prediction (solid line). In the presence of transverse 
anisotropy, it is difficult to determine directly the resonance 
positions from atom loss. Therefore, experimentally the posi- 
tions were obtained by determining the associated atom num- 
ber minima and then subtracting a constant offset. In accord 
with this procedure, we have uniformly offset the experimen- 
tal data in such a way that the measured resonance position in 
the symmetric limit (77 = 1) is equal to the known theoretical 
prediction by Olshanii lO]. The dashed and dot-dashed lines 
are the resonance positions, obtained respectively by solving, 
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/ay 



"irj/C and a\jj /a-y ~ 



See the text for 



more details. 
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